Spin symmetry breaking in bilayer quantum Hall systems 
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Based on the construction of generalized Halperin wave functions, we predict the possible exis- 
tence of a large class of broken spin symmetry states in bilayer quantum Hall structures, generalizing 
the recently suggested canted antiferromgnetic phase to many fractional fillings. We develop the 
appropriate Chern-Simons theory, and establish explicitly that the low-lying neutral excitation is a 
Goldstone mode and that the charged excitations are bimerons with continuously tunable (through 
the canted antiferromagnetic order parameter) electric charge on the individual merons. 



Recently, a canted antiferromagnetic (CAF) state has 
been predicted to exist in bilayer quantum Hall (QH) 
systems at the special filling factor v = 2 M, or more 
generally at v = 2/m where to is an odd integer 
The original theoretical prediction H based on a micro- 
scopic Hartree-Fock calculation has been followed up by a 
number of subsequent theoretical works using a quantum 
nonlinear sigma model || , a bosonic spin approach [hL|J , 
and more detailed Hartree-Fock calculations f^Q] . Fairly 
persuasive experimental support for the CAF phase in 
v = 2 bilayer QH systems also exists |@]. The basic 
idea underlying the novel CAF phase is that the com- 
petition between interlayer tunneling, Zceman splitting, 
intralayer Coulomb interactions, and interlayer Coulomb 
interactions can cause spontaneous symmetry breaking 
in bilayer systems, leading to the CAF phase. This CAF 
phase lies in between the usual spin polarized ferromag- 
netic phase and the symmetric paramagnetic (or singlet) 
phase. 

In spite of the extensive theoretical work on the prob- 
lem using Hartree-Fock or related spin operator ap- 
proaches [|l]-0 , a fundamental understanding of the pre- 
cise nature of the CAF phase, either from the perspec- 
tive of actual QH wave functions or from a long wave- 
length field theoretic viewpoint, is still lacking. In this 
Letter we construct a microscopic wave function for the 
ground state of the CAF state, and develop a Chern- 
Simons theory to study the excitations above the CAF 
ground state. We find a neutral Goldstone mode associ- 
ated with the breaking of the spin symmetry in the CAF 
phase and novel bimeronic charged excitations, which we 
discuss below. Furthermore, we establish that the type 
of symmetry breaking characterizing the CAF state is 
quite generally allowed in bilayer QH systems and may 
in principle exist for a large class of QH states far beyond 
the originally predicted v — 2/m filling factors. 

As pointed out first by Wen and Zee ||, the Halperin 
(to, to, to) wave functions |Io| (neglecting electron spin) 
in bilayer QH systems have the property of fixing the to- 
tal filling factor v in the system, but not the individual 
filling factors of each layer. This allows one to construct 
wave functions that are a superposition of states with 
different numbers of particles in each layer (but fixed 
total number of particles), which leads to spontaneous 



interlayer coherence in the absence of interlayer tunnel- 
ing |B 11 1 . Wc show that the analogous Halperin- type 
construction for spinful electrons in bilayer QH systems 
leads to a CAF phase that breaks the spin symmetry 
spontaneously. 

The electron Hamiltonian for a bilayer QH system can 
be written as 
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Here a, b and a, (3 are layer ("isospin") and spin in- 
dices, respectively; u\ is the intralayer and Uq is the 
interlayer Coulomb interaction; Ag^s is the splitting be- 
tween symmetric and antisymmetric states; A z is the 



Zeeman splitting. In the case when vfj 



(i.e. d = 0, 



where d is the interlayer distance) and As^g = 0, the 
system has a full isospin SU{2) symmetry generated by 
T q = J d 2 x^ a(T T q ab ^ b(J . When only one of d or A S as is 
finite, the isospin symmetry reduces to a smaller U(l) 
symmetry generated by either T z or T x . It is then rea- 
sonable to define single particle wave functions having 
certain transformation properties under the explicit sym- 
metry of the Hamiltonian. For T z symmetry such wave 
functions would describe electrons in a specific layer; for 
T x they are symmetric and antisymmetric combinations 
of electrons in the two layers. In general (and in the rel- 
evant experimental systems) both d and Ag^s are finite. 
Therefore, any U(l) symmetry is destroyed, and there is 
no natural choice of single electron wave functions. We 
choose to work with the symmetric/antisymmetric elec- 
tron wave functions, since in the absence of a gate volt- 
age between the layers, ^ aa-T z b ^ bo does not acquire an 
expectation value and has smaller effect on single elec- 
tron states than interlayer tunneling. Our approach is 
well justified as long as Ag^g is not much smaller than 
Uj — Uq. In the presence of interlayer tunneling, there is 
no degeneracy between different single electron states. Of 
the four possible single particle states, the symmetric spin 
up (S t) state always has the lowest energy and antisym- 
metric spin down (A I) has the highest energy. Since the 
symmetric spin down (5 J.) and the antisymmetric spin 
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up (A |) states may be close in energy, Coulomb inter- 
actions may lead to considerable mixing between them. 

We now construct a Halperin-like wave function for 
our spinful bilayer system. Our wave function does not 
fix the number of electrons in the S j and A | states in- 
dividually, but fixes their sum. Let us first assume that 
the highest energy A j state is completely empty [jl2[ . 
Then, if we label S | states by z, S j states by u, and 
A | states by m, we can easily write the Halperin wave 
function that fixes the number of electrons in the S J, and 
A | states together, but not in each of them separately: 

*({z}{u}{w}) = H( Zl - z 3 r \\izi - Wj y - Uj y 
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Here n and m are odd integers, and I can be any inte- 
ger. A simple calculation then gives the total filling of 
this wave function: v — (n + m — 2l)/(nm — I 2 ). What 
is remarkable about such a wave function is that, since 
the individual filling factors in the S J, and A 1 states 
are not fixed, we can consider wave functions that are 
a superposition of states with various N u — N w . They 
mix states with different values of S z (the z component i 



of spin lies along the direction of the magnetic field) and 
therefore describe states with spontaneously broken spin 
symmetry — the CAF state of It is easy to see 

that taking I = and n = m in (0) gives v — 2/m, i.e. 
the CAF state discussed in In the CAF phase, the 
electrons in the two layers have the same z component 
of spin but opposite x — y components. The direction 
of the Neel order parameter (defined as the difference 
in the spin expectation values in the two layers) comes 
from the spontaneous breaking of the S z spin symmetry. 
It should be mentioned that Halperin wave functions for 
spontaneously broken spin symmetry states may also be 
constructed for single layer QH systems, leading to the 
possibility (at least in principle) of exotic spin states in 
a single layer QH system [ O . 

Properties of the state (g) are conveniently discussed 
using a bosonic Chern-Simons theory [jl4|,[l5|. For sim- 
plicity we again assume that the A J, states are empty 
and consider only three kinds of electrons: for S j, 
* 2 for S |, and * 3 for A f. Eq. (|) tells us that the 
electron is seen as a vortex of strength n by other 
electrons and a vortex of strength I by electrons ^2,3; 
electrons ^2,3 are seen as vortices of strength m and I 
by electrons 1 $>2,3 and respectively. We are therefore 
led to consider the following (bosonic) Chern-Simons La- 
grangian 
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where £cs(a) = j^^ vX diJ,d v ax, and L is a layer index ("Top" or "Bottom") in the Coulomb interaction term. We 
decompose the <IVs into an amplitude, a trivial phase, and a vortex part p5|Jl6[ : = ^fpie %ei 4>vi and ^ a = 
y/p2e ie2 <f) V 2Z a -i for a = 2, 3, with the constraints 4> v i4> v i = <\> V 2 < \>v2 = z a z a = 1. Then (||) can be written as 
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Here Ki — m/pi, terms with "f a b come from the exchange part of the Coulomb interaction, and in the direct part of 
the Coulomb interaction we keep only the layer symmetric part of u\t, which does not vanish in the limit of d = 0. 
By integrating out 6>i and 62 we find that J p = (pi, J) and J M = (p2, J) are conserved. Therefore, we introduce dual 
gauge fields, 6a and b\, such that J M = ^e til/X d v b\ and J M = ^^ vX d v b\. Then, we integrate out the statistical gauge 
fields, a M and a M , and the time component of the dual gauge fields, 60 and 60 • This gives (up to irrelevant constants) 

c = ib a r a + il a (jz + Jfj + {dob a f + ^-i^r (a 6 Q ) 2 + ^- (\Bz\ 2 + {zdzf) + ie al3 b a d b„ + ie^l a d Q l 
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(5) 
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where we have defined vortex and skyrmion currents as 
in @: j; = (J%,JZ) = J- e ^ x d v ^ vl ^ vl ), J* = 

(Jo Ja) = ^ uX du (^2%^ 2 ) , and J* = (j§, J%) = 



±e^d v (i 
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In the ground state there are no vortices or skyrmions, 
so the cancellation of the long range logarithmic interac- 
tion gives two conditions 
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i af3 d a AT + le ap d a b(3 + me^djp] = . (6) 



Recalling that l/2ire a/3 d a bp gives the density of elec- 
trons and l/2ne a Pd a bf3 gives the density of ^ and $3 
electrons, we realize that Eq. (|^) gives us the same filling 
fractions as the Halperin wave function (Q). 

From the last line of (|5|) it is obvious that as we change 
the strength of the Zeeman interaction and/or interlayer 
tunneling, we will stabilize various values of \z\\ and \z%\. 
The most important observation is that when d ^= we 
have 712 > 711 +722, so there is never a direct transition 
from |zi| = 1 & \z 2 \ = to \zi\ = & \z 2 \ = 1, but there 
is always an intermediate phase where both \zi\ = cos#o 
and 1 22 1 = sin 0o are finite. This corresponds to the CAF 
phase for fractional fillings. In this phase interactions fix 
the absolute values of z's but not their relative phase. 
Therefore, when z\z 2 develops an expectation value, we 
have a spontaneous breaking of the U(l) symmetry and 
the appearance of a Goldstone mode. 

In the CAF phase, dynamics of the spin is determined 
by 
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In the CAF phase only (ziz 2 ) develops a non-zero expec- 
tation value, but not (zi) or (z 2 ). Therefore, we can write 
21 = |zi|e l ^+x) and z 2 = \z 2 \e l ^-' l,+ ^ . <j>, the relative 
phase between z\ and Z2, acquires an expectation value 
and gives rise to the Goldstone mode associated with the 
symmetry breaking. We also introduce q — \zi\ 2 — |z 2 | 2 , 
where far away from a vortex core q m in = cos(2#o)- Us- 
ing that fluctuations of q are massive and their gradients 
may be neglected, we find from (0) 
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with a = (711 + 722 — 27i2)/4. Integrating out Sq we get 
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(9) 



We see that the spin wave velocity is v s = (711 + 722 — 
2712) 2 sin 2 (2^o) / (8mp2)- By introducing an infinitesimal 
external Zeeman field and integrating out fluctuations 
in </>, one can also calculate the S z correlation function 
which explicitly shows a Goldstone resonance 
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When the ground state of a system breaks a U(l) sym- 
metry spontaneously, vortices of the U(l) phase will be 
the elementary excitations in the system. For the (1,1,1) 
states discussed in [pp]] , these elementary excitation are 
merons which were shown to have fractional charge [fl"l| . 
However, a pair of merons with opposite vorticity always 
add up to integer charge, or 1. For the wave functions 
defined in (g), one can also imagine a meron in which 
the direction of the Neel vector points outward from the 
center of the vortex, and ask whether such a meron will 
carry an electric charge jl7j . It is clear from the discus- 
sion above that such a meron corresponds to a vortex of 
the z field. Far away from the vortex core, \zi\ = cos#o 
and \z 2 \ — sin 6$, and the relative phase between the two 
z's has nontrivial winding characterized by integer vor- 
ticity, n v . In order to avoid a singularity of this phase, in 
the vortex core we must have either \z\\ = 1 & \z%\ = 
(S vortex) or \zi\ = & |z 2 | = 1 (T vortex). According 
to the definition of Jq , this implies a nontrivial skyrmion 
winding number 



n v x sin #0 S vortex 
— n v x cos 2 #0 T vortex 



fll) 



where n v is an integer characterizing the winding of the 
relative phase between z\ and z 2 . From (^|) the ex- 
tra skyrmion charge has to be compensated by electric 
charge. To cancel the long-range forces, J[—ne al3 d a bp — 
le^dabp] = and J[-le al3 d a bp - me a ^d a l ] = -2ttQ s , 
which immediately gives us the total charge of the meron 
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27 
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(12) 



For 1 = 0, which includes the v = 2/m states discussed 
in 0], we find Q me ron = 1/m x Q s . Note that if we were 
to create simple quasiparticles by squezing a vortex into 
the ground state, Jq = S 2 (x — xq) or Jq = S 2 (x — xq), 
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we could use the same arguments to find their charges: 
q = (m — /)/ (nm — I 2 ) and q = (n — l)(nm — I 2 ). So, as 
in the simple case of a meron in the (m, m, to) state, two 
merons add up to a charge of or the charge of a single 
quasiparticle, q ]l8| |. 

In the simplest case of n = m = 1 and / = in (Q) 
(i.e. ^ = 2), one can give a simple picture of the meron 
excitation in the CAF phase using a generalization of the 
Berry's phase argument in O]. As suggested in H, the 
CAF phase can be described by combining pairs of elec- 
trons into hard core bosons and writing the wave function 
as |*) = cos6\S) +e**sm6\T). Here IS") and \T) denote 
singlet and triplet bosons, respectively; the relative phase 
between the two bosons, <f>, determines the direction of 
the Neel vector in the x — y plane. When a vortex is 
present, this phase winds nontrivially around the vortex 
core and is characterized by an integer vorticity, n v . At 
the center of the core, one has to demand that there is 
only one kind of boson present (so as to avoid a singu- 
larity of the relative phase); therefore, one expects the 
appearance of two kinds of vortices: vortices with a sin- 
glet core (S) or a triplet core (T). We can now imagine 
taking a pair of electrons and adiabatically moving them 
around the vortex. In the course of such adiabatic trans- 
port, the wave function for a pair of electrons will acquire 
a phase iT — §{ijj\dvb) = i n„ sin 2 for an S vortex or 
iT = —i n v cos 2 for a T vortex. The Berry's phase in 
adiabatic transport is indistinguishable from extra flux 
going through the system A$ = <&o/2 x r/(27r), where 
the factor of 1/2 comes from the fact that we transported 
a pair of electrons. This extra flux can be related to the 
charge carried by the meron as Aq = a xy A$ = n v x sin 2 
and Aq = —n v x cos 2 9 for S and T vortices, respectively. 
So, the two kinds of merons in this case carry fractional 
charge; the charge depends on where the system is in the 
phase diagram, i.e. on the CAF phase order parameter 
9 ( 9 goes to at the boundary of the CAF phase with 
the spin singlet phase, and n/2 at the boundary with the 
fully polarized ferromagnetic state; see || for details ). 
However, two merons with opposite vorticities again add 
up to a charge of or 1, as in the v = 1 bilayer (1, 1, 1) 
state 0. 

It is also instructive to consider an explicit wave func- 
tion for a meron in the CAF phase at v = 2. As 
discussed above, the wave function of the CAF phase 
may be conveniently written as (in the limit when d is 
small) |*o) = I!™ ( c °s Sj, + e^sin0 T^)|0), where 

&m = l/v / 2( c SmT C Sml _ C Sml C SmT) an °- = C Smt C Am1 

create singlet and triplet combinations of electrons with 
orbital momentum to, and |0) is the Fock vacuum. Us- 
ing the definition of the Neel order parameter, N a (z) = 

(*|S?-5S|*)=E TO n*^(*)*n(*)x(*l4ma^/9^- 

c Ama< T af3 c Sni3\^) , where * m (z) is the wavefunction of an 
electron in the first Landau level with angular momen- 
tum to, one can easily prove that state |$o) has a uniform 



N in the XY plane N+(z) = 1/2 cos sin 0e^. To have a 
meron we need a wave function where the direction of the 
Neel vector winds around as one goes around the center of 
the meron. This is achieved by considering the following 
wave function |* M ) = U m ( cos9 S l + e^sinfl ft,)|0), 
where f} n = c Sm+1 ^c' Am ^. For |*m) one finds that 

N+(z) = l/2E m *mri*m+i (z)cos0sin0e^°. Since 
tx z m exp(— \z\ 2 ) we find that a.rcta.n(N y / N x ) = 
arg(z) + (j>; the direction of N winds in the XY plane 
following the argument of the complex coordinate z. It 
is also obvious from |^m) that it describes a state with 
a missing electron in the S f state of m = 0, so we have 
an S vortex with charge — sin 2 0. 

In summary, we have developed an analytic theory for 
the bilayer QH CAF phase. Our theory is consistent with 
the original Hartree-Fock theory for v = 2, but is gen- 
eral enough to predict a whole new class of fractional QH 
CAF phases as well as the correct excitation spectra. 
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